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We study the ground-state phase diagram of a Bose-Fermi mixture loaded in a one-dimensional optical lattice 
by computing the ground-state fidelity and quantum entanglement. We find that the fidelity is able to signal 
quantum phase transitions between the Luttinger liquid phase, the density-wave phase, and the phase separation 
state of the system; and the concurrence can be used to signal the transition between the density-wave phase and 
the Ising phase. 



Ultra-cold atomic gas loaded in optical lattice are attracting 
more and more attentions due to ambitions of getting deep 
insight into some essential physical phenomenons, such as 
quantum phase transitions (QPTs), 1 in the condensed matter 
physics. The prediction 2 and the successful observation^ of 
the QPT from a superfiuid to a Mott-insulator was one of the 
greatest works in this field. Quite recently, experimental pro- 
gresses are very promising for studying more non-trivial quan- 
tum phases in clod atomic systems. For example, bosonic 
and fermionic atoms can simultaneously be trapped in opti- 
cal lattice in a controllable way. This ultra-cold atomic sys- 
tem, the so-called Bose-Fermi mixture, often reminds peo- 
ple of the solid state systems including the electron-phonon 
interaction. The latter have been studied for a long period 
and have a complicated quantum phase diagram. So to ex- 
plore new possible phases in this system becomes an interest- 
ing theoretical problem. Along this line, several works had 
been done. 4,5 ' 6 ' 7 ' 8 It is worthwhile to mention that Lode Pollet 
et. al. recently studied the ground-state phase diagram of a 
Bose-Fermi mixture loaded in a one-dimensional (ID) optical 
lattice by using quantum Monte Carlo (QMC) simulations. 8 
Several phases, including Luttinger liquid (LL) phase, den- 
sity wave (DW) phase, phase separation (PS) state, and Ising 
phase, are predicted. 

In recent years, some concepts emerging in the quan- 
tum information theory 9 are extensively used to study the 
critical phenomena in quantum many-body systems. One 
of the typical examples is the entanglement. Many efforts 
have been made to the relation between the entanglement 
and QPTs. 101112 Quite recently, the fidelity, as a measure 
of similarity between states, was proposed to study the crit- 
ical phenomena. The motivation is very simple: a dra- 
matic change in the structure of the ground state around 
the quantum critical point should result in a great differ- 
ence between the two ground states on the both sides of 
the critical point. The fidelity has been successfully ap- 
plied to study the spin, fermionic, and most recently bosonic 
systems i 13 ' 14 ' 15 i 16 ' 17 ' 18 ' 19 ' 20 ' 21 Compared with entanglement, 
the fidelity is purely a geometrical quantity; an obvious ad- 
vantage is that in analyzing the QPTs it does not require a 
priori knowledge of the order parameter and the symmetry of 
the system. 

In this paper we try to study the ground-state fidelity and 



the entanglement of the Bose-Fermi mixture in a ID optical 
lattice. The aim is two-folded: one is to test the role of fi- 
delity and entanglement in a more realistic system, the other 
is to study the ground-state properties of the Bose-Fermi mix- 
ture. Follow Lode Pollet et. al}, we assume that a mixture 
of bosonic and fermionic atoms is loaded into a ID optical 
lattice, and the temperature is low enough such that quantum 
degeneracy is achieved. The system is then described by a 
lowest-band Bose-Fermi Hubbard model, 
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where b, (bf) and c, (ct) are the bosonic and fermionic an- 
nihilation (creation) operators at site i, respectively. Bosons 
(fermions) can hop from site i to the nearest neighbor site i + 1 
with tunneling amplitude ?b (?f). Furthermore, a large occu- 
pation of bosons on a single site is suppressed by the on-site 
repulsion interaction Ubb- Bosons and fermions can mutually 
repel or attract each other on each site depending on the sign 
of Ubf- In this paper, we choose Ubf > 0, and consider the 
case where both of the bosons and the fermions have a density 
as: Nf = Nb = N/2. 

We now briefly introduce the ground-state phase diagram 
of the model. When the Ubf is small enough, the fermions 
behave as a LL, the interaction between them is induced by 
the bosons. At the same time, the bosons form an interacting 
liquid too. So we have a LL of fermions, which weakly inter- 
acts a boson liquid. When Ubb is small and Ubf is large, the 
system is in the first order unstable to the PS with hard domain 
walls; in other words the system is separated into two regions: 
a bosonic region and a fermionic region. When both Ubb 
and Ubf are very large, the bosons behave as fermions, which 
means that the occupation of more than one boson at a single 
site is not allowed. The model can then be mapped into the ID 
XXZ model: H X xz = £/ ^(ofcrf+i +o y i o y i+1 )+J z a^a z i+v where 
J = -(?Bf F )/t/BF and P = (f 2 + f 2 )/(2t/ BF ) - f 2 /(2t/ BB )^. 
This implies that there are three phases in this limit: the fer- 
romagnetic phase, a gapless DW phase and a gapped Ising 
phase. In the ferromagnetic phase, boson-boson bonds and 
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FIG. 1: The obtained fidelity F(26Ubf, Ubf, Ubb) when tp = 
4.0, fa = 1.0, SUbf = 0.5. It can be observed obviously that there 
are two phase transition boundary lines indicated by the drop in the 
fidelity. The phases such as LL, DW, and PS are identified by the 
comparing this phase diagram with the one which was proposed ear- 
lier by the calculation of the correlation functions (see text). N = 8 
(APBC). 

fermion-fermion bonds are favored compared with bonson- 
fermion bonds due to the larger exchange interactions, this 
mechanism makes the system form two regions with hard do- 
main walls. So the ferromagnetic phase corresponds to the 
PS in the mixture. This mechanism is similar to the one ap- 
pearing in the ID asymmetric Hubbard model 23,24 where the 
PS is also occurred when the system is away from half-filling. 
While in the DW phase and Ising phase, the system always 
favors bonson-fermion bonds. We could like to emphasis that 
the PS in the large Ubb limit and the one in the small Ubb 
limit are different since that the latter allows the occupation of 
more than one boson on a single site. In the whole PS region, 
with the increasing of Ubb, the boson repulsion exerts a pres- 
sure such that the region occupied by the bosons will grow 
and at the same time the local density of bosons will decrease. 

As mentioned, the fidelity is nothing but the modulus of the 
overlap of two ground states relative to two different choices 
of the Hamiltonian parameters. In this paper, we mainly focus 
on these two: 

F{25Ubf, Ubf, Ubb) - KiA(7 B F-*f/BFl | Af/BF+*f/BF)l' (2) 
F(26Ubb, Ubf, Ubb) - K'Ac/BB-tff/BBl'Ac/BB+tft/BB >l, (3) 

in which \iffj) stands for the ground state of the Hamiltonian 
(1) with the parameter A, and is calculated by the Lanczos 
method for a finite sample. To avoid the ground-state level 
crossing, anti-periodic boundary conditions (APBCs) are ap- 
plied for system size N = An and periodic boundary condi- 
tions (PBCs) for N - 4n + 2, where n is an integer. According 
to the original motivation of the fidelity, a drop in the fidelity 
of two ground states separated by two lightly different param- 
eters is expected to be a signature of the QPT. 

In Fig. 1, we show one of our main results, i.e. the fidelity 
F(25Ubf, Ubf, Ubb) defined on the Ubb - Ubf plane. Com- 
pare this figure with Fig. 4 of Ref. 8, perfect similarity can 
be observed. The boundary lines of phase transitions between 
the LL, the DW and the PS are clearly indicated by the drop of 



the fidelity. We would like to emphasis that the phase diagram 
presented here is obtained in such a small cluster and with- 
out any knowledge of the correlation properties of the system. 
It is also clearly observed that the drop of the fidelity along 
the phase transition line between the DW and PS phases be- 
comes deeper and deeper as the interaction decreases. This 
phenomenon indicates that although the phase transition is 
within the same class but the similarity of the ground state 
is changing along this line. 

However, the phase transition, as reported in Ref. 8, be- 
tween the DW and the Ising phases is not indicated in Fig. 1. 
According to the effective model, i.e. ID XXZ model, this 
transition belongs to the KT universality class. 8,25 Lode Pollet 
et. a/. 8 did numerical calculation of the correlation functions 
in a relatively larger system (N ~ 30) by the QMC simula- 
tions and claimed that a true long-range order may exist in 
the Ising phase. It is highly difficult for us to make a scaling 
analysis of the correlation functions by the Lanczos method. 
However, the ID XXZ 26 model provides us a clue to investi- 
gate this problem. It was reported that the concurrence,— as a 
measure of entanglement between two qubits, reaches maxi- 
mum at the SU(2) point 11 of the XXZ model. This maximum 
point corresponds to the transition point between the DW and 
the Ising phases. 

The concurrence in the spin models can be calculated in 
the following way. Due to the global SU(2) symmetry of the 
XXZ model, the z-component of the total spin of the system 
is a good quantum number, the reduced density matrix p J]!+ i 
of two neighboring spins has the form 

( u + \ 



{ u ) 

in the spin basis | TT), I IT), I ti), I ID- The elements in the 
reduced density matrix py+i can be obtained from the corre- 
lation functions 

u ± =^(\±2{o*) + {cj z i o z M )), (5) 
w= (6) 

z = \((cr?<T? +1 ) + + i<ofo? +1 > - /<«,»• (7) 

All the information needed is contained in the reduced density 
matrix, from which the concurrence is readily obtained as ? 

C = 2 max [0, \z\ - Vm + m ] , (8) 

which can be expressed in terms of the correlation functions 
and the magnetization by using equation (5), (6) and (7). 

The concurrence is only valid for two-qubit system. For 
the Bose-Fermi Hubbard model, the double occupation of two 
particles is almost not allowed, the state of the two neighbor- 
ing sites can be described by four basis \ff), \bf), \fb), \bb) 
where / (b) represents there is only one fermion (boson) on a 
single site. So we can associate a pseudo-spin "up" ("down") 
with the "/" ("&")• Then the equations described above can 
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FIG. 2: The obtained concurrence in the DW phase where fp = 
4.0, ?b = 1.0. The concurrence of the N = 6 (N = 8) had been 
reduced by 0.029 (0.009) in order to make the results reside in the 
same region in figure. All other parameters are indicated in the fig- 
ure. 



also be perfectly applied to the Bose-Fermi Hubbard model. 
In other words, if we calculate the trace of the reduced matrix 
of two nearest sites, which reads T - Trp;,i + i, then T must 
equal to 1 in this limit. While T only approximately equals 
1, when Ubb and Ubf are not infinite but large. As long as 
(1 — T) is small enough the concurrence calculated by Eq. (8) 
is a good characterization of entanglement between two sites. 
In this way, we are able to calculate the concurrence in the 
Bose-Fermi mixtures approximately. This kind of treatment 
was also used by some other group recently 28 . Obviously, as 
expected a maximum is clearly observed in the behavior of 
the concurrent, which indicates a QPT point. It is important 
to point out that the variety of the concurrence is not a dra- 
matic one in the whole DW region. Using the exact solution 
of the XXZ chain, the phase transition from the DW phase to 
the Ising phase should occur at —J = J z , i. e. Ubb = 6.7 
for f F = 4.0, f B = 1 -0, Ubf = 60. But according to the ob- 
tained concurrence, the max point U^ x approximately equals 
to 7.7. The discrepancy comes from the approximations made 
and the size effect. 

Some skeptical readers may wonder why other kinds of 
measurements of the entanglement, for example the von Neu- 
mann entropy, are not calculated here since these kinds of 
measurements are more suitable for the Bose-Fermi Hubbard 
model at the first glance of their definitions. We actually cal- 
culated von Neumann entrpy though we did not show here. 
But it turns out that these quantities are more difficult to wit- 
ness the phase transition between the DW and Ising phase. 

After all the analysis, we contribute the missing of the phase 
transition signature between the DW phase and the Ising phase 
in the fidelity to two reasons. The first one is that this kind 
of transition is actually a very weak one, which means that 
the change of the ground state around the critical point is not 
dramatic, at least in the finite size system according to our 
results. The second one is that, as reported before, the fidelity 
may not be a good indicator of those transitions of infinite 
order, such as KT transitions. 17 However, further studies are 
definitely needed in order to answer the question completely! 

We also calculated the fidelity F (2SUbb, Ubf, Ubb) which 
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FIG. 3: The obtained fidelity F(2SU SB , Ubf, Ubb) (a), (b) and the 
first order derivative of the local density of boson (see text)(c), (d) 
when tp = 4.0, ?b = 1.0, Ubf = 60.0. All other parameters are indi- 
cated in the figure. 



can signal the transition between the PS phase and DW phase 
too. As you can see in the Fig. 3, the most dramatic drop 
in the fidelity is in correspondence with the phase transition 
point between the PS state and the DW phase. The critical 
point find here is consistence with the one found in the Fig. 1 . 
In addition, some other drops in the fidelity can be observed 
in the PS region. These drops are related to the changing rate 
of the local density of boson. To show this, we define the local 
density of boson as, 



1 N 



1)>, 



(9) 



and its first order derivative is calculated. Comparing Fig. 3 
(a) and Fig. 3 (b) with Fig. 3 (c) and Fig. 3 (d) respectively, 
every drop in the fidelity has its company, a drop in the deriva- 
tive of the local density of boson. This relation can be under- 
stood by the dramatic change of the local density of boson 
will lead to a big change of the ground state wave-function, 
and then will make the fidelity decrease greatly. The drop of 
the fidelity at the QPT point between the PS state and DW 
phase can also be thought like this way. These observations 
strongly imply that the transition between PS phase and DW 
phase is within the Landau's symmetry breaking theory and 
a first order one. Furthermore, it is easy to notice that the 
phase transition point is not affected by the system size and 
the boundary conditions used, but the drops found in the PS 
phase is strongly affected, which indicates that these drops 
may be a size effect and can not be identified as phase tran- 
sitions. Furthermore, the transition between DW phase and 
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Ising phase is not observed again! 

In summary, we calculated the fidelity F(26Ubf, Ubf, Ubb) 
and F(26Ubb, Ubf, Ubb) of the ID Bose-Fermi Hubbard 
model, which can be used to describe low-temperature physics 
of the atomic Bose-Fermi mixtures loaded in ID optical lat- 
tices. It is showed that although the ground state phase dia- 
gram of this system is complicated, and the fidelity may be a 
good tool to study the phase transitions without a priori knowl- 
edge of the order parameter and the symmetry of the system. 
But one should be very careful because the fidelity may fail for 
the transition which is belong to the KT-like type, for exam- 
ple, the transition between DW phase and Ising phase in this 
system. We also calculate the concurrence in the DW phase. 



The result indicate a QPT may exist in the DW phase although 
the fidelity have no singulary at this transition point. 
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